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ABSTRACT 

We construct a sample of 75,863 star forming galaxies with robust metallicity and star forma- 
tion rate measurements from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), from 
which we select a clean sample of compact group (CG) galaxies. The CGs are defined to be 
close configurations of at least 4 galaxies that are otherwise apparently isolated. Our selection 
results in a sample of 112 spectroscopically identified compact group galaxies, which can be 
further divided into groups that are either embedded within a larger structure, such as a cluster 
or large group, or truly isolated systems. The compact groups then serve as a probe into the 
influence of large scale environment on a galaxy's evolution, while keeping the local density 
fixed at high values. We find that the star formation rates (SFRs) of star forming galaxies in 
compact groups are significantly diff'erent between isolated and embedded systems. Galaxies 
in isolated systems show significantly enhanced SFR, relative to a control sample matched in 
mass and redshift, a trend not seen in the embedded systems. Galaxies in isolated systems ex- 
hibit a median SFR enhancement at fixed stellar mass of -h0.07 ±0.03 dex. These dependences 
on large scale environment are small in magnitude relative to the apparent influence of local 
scale efl'ects found in previous studies, but the significance of the diff'erence in SFRs between 
our two samples constrains the eff'ect of large scale environment to be non-zero. We find no 
significant change in the gas-phase interstellar metallicity for either the isolated or embedded 
compact group sample relative to their controls. However, simulated samples that include ar- 
tificial off'sets indicate that we are only sensitive to metallicity changes of log 0/H >0.13 dex 
(at 99% confidence), which is considerably larger than the typical metallicity differences seen 
in previous environmental studies. 

Key words: galaxies: interactions, galaxies: groups: general, galaxies: evolution, galaxies: 
abundances 



1 INTRODUCTION 

The evolution of a galaxy is fundamentally shaped by its envi- 
ronment. Qualitative evidence of the influence of environment has 
been known since the 1930s (e.g., [Hubble & Humason|T93Tl ); the 
large surveys now available have added their statistical weight to 
its influence, allowing a quantification of environmental changes. 
Nearly every observable quantity of a galaxy has been shown to 
vary with environment, including star formation rates (SFRs), mor- 
phologies, colours, active galactic nucleus (AGN) fractions, and 
mean stellar mass (e.g., Dressier" 1980', 'Balogh et al.^ 2004||Kauff-| 
mann et al. 2004 , Croton et al. 2005 , Baldry et al. 2006[ |Park et" 
2007i r The primary mechanisms that can operate in denser envi- 
ronments are strangulation, harassment, and ram-pressure stripping 
(e.g., |Gunn & Gott|1972[ [Larson et al.|1980[ jMoore et al.p996l ). 
However, disentangling which of these might be a fundamental 
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variable and which are corollary effects is a non-trivial task (e.g., 
iWhitmore & Gilmore.1991 , Blanton et al.[2005[[Blanton & Berlind[ 
[2007 1 [Skibba et"aL]|2009j ). Moreover, relatively simple processes 
such as individual galaxy-galaxy interactions become complicated 
by the influence of a local, high density environment. For example, 
at low densities, close galaxy pairs show enhanced SFRs ( Lam-[ 
[bas et al.|200"3l|Alonso et al.|2004[rEmson et al. 2010) . However, 
once the galaxy pair is found within a high density environment, 
the strength of the SFR enhancement diminishes or becomes unde- 
tectable ( jAlonso et al.l2004 2006, B aldry et al. 12006 ) [E lHson etaL] 
[2010[ ). Clearly, galaxies are sensitive to their environment beyond 
the scale of their nearest neighbour. The question remains, how far 
beyond the nearest neighbour does environment matter? Is it merely 
a local density effect, or does cluster membership, with scales far 
beyond the local, also impact a galaxy's evolution? This work will 
attempt to address this question. 

Several density trends are well documented in the literature for 
low redshift galaxy samples. The morphology-density and colour- 
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density relations describe the declining fractions of blue, late-type 
galaxies with increasing density 0ressler 1980, Postman & Geller 



1984 



2004 



Balogh et al.|1997| |1998|[Martmez et al..2002, Tanaka et al. 
•61 



Baldry et al. 2006)1. The overall distribution of star formation 



rates also correlates with density, with the average galactic SFR 
declining a s density increases (Hashimoto et al. 1998 , Poggianti 
let al.|2008[[Kauffmann et al.|2o64[ [Cooper et al.|2008| ). Interpreta- 
tions of the SFR-density relation are complicated by the increasing 
fraction of low- (or non-) star forming elliptical galaxies as density 
increases. The contribution of early type galaxies will bias the av- 
erage SFR to lower values at higher densities (e.g., |Balogh et al.| 
[2004] ). 

The dependence of SFR on density becomes much less clear 
for only the star forming subsample of galaxies, with many conflict- 
ing interpretations. Some studies find that among the star-forming 



fraction, there is no residual dependence on density ( Cou ch et al. 
200n [Balogh et al.||2004[ [Tanaka et al.|[2004 Weinman n et al. 



20W 



Patiri et al."20b6 . Park et al.|2007|[Peng et al.|2010||McGee 



et al. 2011, Ideue et al. 2011}. According to these works, the ob 



served decrease in SFR with increasing density is entirely due to 
the changing fraction of star-forming galaxies with density, and 
once this dependence is eliminated, the distribution of SFRs is in- 
dependent of density. Other studies counter these results, finding 
that the population change is not enough to account for the en- 
tire SFR-density relation, and that a density dependence remains 

flogh et al.|1998[ [Pimbblet et al. 12002] [Gomez et all2003l|We^ 
la et al. 2008 ). Moreover, the scale over which the SFR is influ- 
enced by its surroundings is widely debated. Many authors find that 
it is primarily on the local scale (generally considered to be of order 
1 Mpc or less) that the residual SFR-density dependence is found 
([Hashimoto et al." 199 8 " Carter et al. "2001", Lewis et al.7002VKauff^-" 
[man n et al. 2004, Blanton et al. 2006, Blanton & Berlind 2007 ). 
By contrast, some studies indicate that the environment on scales 
of order several Mpc is of greater importance than that of sub-Mpc 
scales ( [Goto et al.|2003l[Park & Hwang|2009[ ). 

A further probe into the scales involved in environmental ef- 
fects comes from the gas phase metallicity, which, like the SFR, 
can be modulated by the gas supply (Ellison et al. 2008 a( [Man-[ 
[nucci et al.|[2010[ [Lara-Lopez et al.[[2010, ,Yates et al.„2011[ ). In 
samples of interacting pairs, changes in the gas phase metallicity 
and SFR are linked, and have been shown to respond to some of 
the same physical processes dSc udder et al.[[2012 ). Theory off'ers 
the picture of an interaction triggering large scale gas flow from the 
outer regions of a galaxy to the centre, and the enhanced densities 
of gas sparking star formation (Mihos & Hernquist 1996, Barnes 
& Hernquist 1996 , Barnes 2004, Rupk e et al.|201 0, Montuori et al. 
2010^ Torrey et al. 2012) . In agreement with this theoretical pic- 



ture, close galaxy pairs show enhanced SFRs (Larson & Tinsley 



fT978',"Donzem & Pastoriza' 1997 , "Barton et al.'2000i"Lambas et al. 
2003 , Alonso et al. 2004 2006 , Woods & Geller 2007 , EUison et al. 
2008b 2010 ) and diluted metaUicities ( Kew ley et al.|2006HEllison[ 
jet al.„2008b, Michel-Dansac et al..,2008j . This indicates that gas 
phase metaUicities ought to change, as the SFRs do, with density. 
In confirmation of this idea, studies of galaxies in clusters and other 
dense environments have shown metal enhancement relative to the 
field ( [Mouhcine et al. [2007 [[Cooper et al.|2008l[Emson et al.[2009[ ). 

In this paper, we employ a new tactic to investigate the eff'ects 
of large scale environmental dependences. We select two samples 
of galaxies with similarly high small scale densities, but with dif- 
ferent large scale environments. To this end, we use the catalogue 
of 3491 compact group (CG) galaxies in 828 CGs presented in 
[Mendel et al.[ ( |20TT[ henceforth Mil) as a refinement on the sam- 



ple of McConnachie et al.[ ( [2009 1 henceforth M09). Although the 
original CG selection criteria were designed to identify groups that 
are not part of a larger overdensity ( jHickson,19 82 , M09), Mil find 
that the SDSS CG sample can be divided into two distinct popula- 
tions. One population appears to be truly isolated, whereas the other 
population appears embedded within a large scale structure, but iso- 
lated from other galaxies within the group. Mil show that despite 
the strong diff'erences in their large scale environments, the galaxies 
in these 'embedded' and 'isolated' groups have similar morpholog- 
ical properties. In this paper, we investigate the spectroscopic prop- 
erties of the embedded and isolated CG galaxies, focussing on their 
star formation rates and metaUicities. The photometric catalogue of 
CG galaxies is therefore cross-matched with a spectroscopic sam- 
ple of emission line galaxies from the SDSS DR7 ( A bazajian et al.[ 
[2009 1 ). In Section |2] we describe the sample selection, including 
metallicity and AGN calibrations. In Section [3] we describe our 
methods for quantifying differences between the isolated and em- 
bedded CG galaxies and a sample of control galaxies. In Section]?] 
we discuss the implications of our results, along with a comparison 
to previous works, and present our conclusions in Section[5] 

We assume Qm = 0.3, Qa = 0.7 and Hq = 70 km s"^ Mpc'^ 



2 SAMPLE SELECTION 

Our spectroscopic sample is taken from the publicly available 
MPA-JHU SDSS DR7 catalogu£]of 927,552 galaxies. This cat- 
alogue provides measurements of up to 12 emission lines per 
galaxy, corrected for stellar absorption lines and Galactic redden- 
ing. Of these lines, we will eventually require robust flux measure- 
ments to be present in Hof, I^g, [Oii]^3727, [Oiii]^^4959, 5007, and 
[Nii]/16584 for metallicity calibration^ As some of these galaxies 
are not unique objects within the catalogue, we use the MPA-JHU' s 
duplicate catalogu^to remove duplicate galaxies from the sample. 
In this section, we describe the quality control measures taken in 
order to select a reliable sample of star-forming galaxies with mea- 
surable metaUicities. 



2.1 Quality Control 

In order to ensure that the flux measurements will be robust for 
subsequent metallicity calculations, careful checks were made on 
the quality of the flux measurements. As an initial step, all galaxies 
with zero and negative flux values were excluded. Flux values of 
zero were determined to be either due to clipped lines or due to the 
line being redshifted out of the spectral range of the SDSS. Neg- 
ative flux values, by contrast, seemed to be largely due to poorly 
subtracted B aimer absorption. 

A redshift cut was imposed on the galaxies to ensure that the 
spectral lines needed for the metallicity calibrations will not have 
been shifted out of the spectral range or into a region of the spec- 
trum dominated by sky line residuals. The bluest wavelength of 
the emission lines needed in our analysis, [Oii]/13727, dictates the 



^ Available at: 'h ttp : //www . mpa- garching . mpg . de/SDSS/DR7/raw_[ 
data.html 

^ Recently, Groves et al. (2011 » have found that the I^ equivalent widths 
from this catalogue are systematically underestimated by 0.35A. We do not 
expect this to affect our results, as all our comparisons are relative, and the 

underestimation fou nd is constant across the entire sample. 

^ Available at ht tp : //www . mpa- garching . mpg . de/SDSS/DR7/Data/[ 
^all_ma tches_dr7 . dat[ 
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lower redshift limit. This line has to be redshifted at least to the 
lower limit of the SDSS spectra at 3800 ^ At the other end of the 
spectrum, significant sky line residuals appear at ~ 8000 A. We use 
Hof to set the upper redshift, by requiring that Hof cannot be shifted 
past 8000 A. These wavelength constraints result in a redshift range 
of 0.02 < z< 0.25. Additionally, we require that the catalogue flag 
z_warning be set to 0; a non-zero value for this parameter indicates 
a bad redshift. 

An inspection of the distribution of flux values for emission 
lines showed that both the flux values and their errors were bi- 
modal with a prominent secondary peak at flux values > 10~^^ erg 
s~^cm~^. A number of the galaxies with the highest flux values were 
visually inspected and the emission line fluxes were re-measured 
manually. The manually measured emission line fluxes were found 
to be within the bounds of the lower flux peak, indicating that the 
high flux values were spurious. The continuum fluxes and errors 
were also found to be bimodal. The secondary peak in the contin- 
uum error (> 10~^ erg s~^cm~^) pushes the continuum fluxes, and 
thus the line flux, to extremely large values. Imposing a cut around 
the main distribution of continuum flux errors eliminated the ma- 
jority of the secondary (high value) flux peak. The remainder of the 
high fluxes were found to be associated with high flux errors, and 
imposing a flux error threshold eliminated the last remnants of the 
secondary peak in the flux values. 

It is standard practice to impose a Signal-to-Noise (S/N) cut 
on emission line fluxes to ensure high quality metallicity determi- 
nations (e.g. Kewley & Ellison 2008). We adopt a fairly standard 
S/N cut of 5 on all emission lines needed for our metallicity cal- 
ibrations. We also impose a S/N> 5 cut on the Balmer ratio (the 
ratio of Hof to H/? flux) to minimize the scatter away from the the- 
oretical Case B recombination limit. Our placement of the S/N cut 
in the Balmer ratio is motivated by the data; less stringent S/N cuts 
yield a rapid increase in the scatter in the Balmer ratio. However, 
increasing the severity of this cut would not significantly decrease 
the scatter, and would serve only to limit our statistics. As we ex- 
pect some random scatter below the theoretical lower limit of 2.85, 
we permit 3cr scatter below the theoretical limit on the Balmer ratio 
(see Figure [T]). This ensures that all Balmer ratios are reasonable, 
and does not exclude a large number of galaxies. A summary of the 
final quality control cuts is presented in Table [T] 

With robust flux values in hand, we next correct for inter- 
nal galactic reddening to obtain intrinsic fluxes, using the Small 
Magellanic Cloud (SMC) extinction curve presented in |Pei|(T992| . 
E(B-V) values for the SMC and the Milky Way are very similar at 
these wavelengths, so our decision to use the SMC curve does not 
aff'ect our results. 



2.2 AGN Removal 

Before metallicities can be calculated, any galaxies containing an 
Active Galactic Nucleus (AGN) must be identified and excluded. 
AGN are known to alter emission line strengths, which makes 
any metallicity calculated for a galaxy hosting an AGN unreliable. 
The canonical method of discriminating between star-forming and 
AGN galaxies is through the use of emission line ratios (Baldwin, 
|et al.||1981[ [Veilleux & Osterbrock]|1987| ). In particular, the most 
broadly used diagnostic plot is the so-called BPT diagram, after the 
[Baldwin et al.,(1981j work in which it was proposed, which plots 




3(7 scatter 

H«/H^ = 2.85 



3 6 9 12 15 

Balmer Ratio (Ha/H/3) 

Figure 1. Final sample Balmer ratios vs. Signal to Noise for the quality 
controlled sample. Blue dashed line is the 3cr lower limit to the theoretical 
limit of 2.85 (red vertical line). 



Final Quality Control 



Redshift 
Continuum error 
Flux error 

Hof/I^ Signal/Noise 
Flux S/N 



0.02 < z < 0.25 
10-19 < < 10^-15 

lO-l'^ < Ferr < lO'^^ 

S/N>5 
S/N>5 



Table 1. Final quality control parameters. Continuum and flux errors are 
in of units erg s-^cm-^. Ha, I^, [Oii]i3727, [Oni]ii4959, 5007, and 
[Nn]/16584 must pass these cuts to be considered a reliable detection in 
that line. 



[Oiii]i5007/Hyg against [Nii]i6584/HQf. The [Oiii] line is colUsion- 
ally excited in both star-forming and AGN regions, so it cannot 
distinguish between the two ionisation sources, but rather serves 
as a metric for the strength of the ionisation field. [Nii], by con- 
trast, originates in partially ionised regions, which are found pref- 
erentially around AGN sources due to the power-law radiation they 
produce. Photoionisation from O & B type stars does not produce a 
large partially ionised region, so the [Nii] line strength is relatively 
small in galaxies dominated by young stellar light. The combina- 
tion of these two line ratios in the BPT diagram thus allows for 
diff'erentiation between galaxies dominated by star formation and 
those with fluxes dominated by an AGN. In order for a galaxy in 
our sample to pass this diagnostic, it must have reliable detections 
in all four of the required emission lines. 

There are a number of possible lines drawn on the BPT dia- 
gram to distinguish between star-forming and AGN. We select the 
diagnostic of |KaufiTnann et al!] ( |2003 b K03). K03 is an empirical 
model, designed to follow the left-hand wing of the BPT diagram 
and is defined as follows: 



log([OIII]^5007/H/?) = 



0.61 



log([NII]^6584/HQf) - 0.05 



+ 1.30. (1) 



4 We note that relaxing the redshift criteria from 0.02 to 0.01 does not in- The K03 demarcation follows the star forming wing closely, 

crease the sample size of our final CG sample. (See Figure 3 of Ml 1.) so the fraction of AGN galaxies included in the SF sample should 
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be low; [Stasinska et aL] ( |2006| ) estimates only 3% contamination 
from AGN galaxies. We classify the sample of galaxies as star- 
forming or AGN based upon this diagnostic, resulting in a final star 
forming sample of 126,756 galaxies. 



2.3 Metallicity Calibrations 

The metallicity calibration adopted here is the adaptation of the 
[Kewley & Dopita,^2002 henceforth KD02) recommended method 
presented in ?, KE08, which we refer to as the KD02-KE08 
method. This calibration does not require the detection of as many 
strong lines (and in particular does not require the [Sii] line), which 
helps to keep the final sample as large as possible. The KD02-KE08 
calibration has low intrinsic scatter, and is easily converted into 
other metallicity calibrations without significant residuals. 

The original KD02 calibration is based on photoionisation 
models and stellar population synthesis for high metallicities, and 
an average of R23 methods for low metallicities, where R23 is de- 
fined as 



^23 — 



[OII]i3727 + [OIII]^4959 + [OIII]i5007 



(2) 



The modification implemented by Kewley & Ellison (2008) 
is in using a diff'erent set of R23 methods for the low metallicity 
galaxies. For high metalhcities (log([Nii]/[Oii]) > -1.2), the KD02- 
KE08 calibration is the same as presented in KD02. This calibration 
requires only [Nii] and [On] detections; as long as the galaxy is 
in the upper branch regime, we do not require any further strong 
line detections for the metallicity calibration. For log([Nii]/[Oii]) 
< -1.2, the lower branch has changed from the original average of 
the |Zarits ky et al. (1994 ) and McGaugh"('1991') calibrations to an 
average of the Kobulnicky & Kewley ( 2004 ) and McGaughj ( |1991| ) 
lower branch calibrations, all of which are based on the /?23 and O32 
diagnostics, where O32 is defined as 



O32 = 



[0III]/14959 + [OII1]^5007 



(3) 



[OII]^3727 

The [Kobulnicky & Kewley | ( [2004| ) calibration takes an iterative ap- 
proach, solving both for the ionisation parameter (defined as the 
number of ionising photons per number of hydrogen atoms) and 
the metallicity simultaneously. As these lower branch calibrations 
require the [Oiii] lines, galaxies that fall in the lower branch have 
the additional criteria that they must have significant detections in 
[Oiii]ii4959,5007. 



2.3.1 Final Metallicity sample 

Our final metallicity sample is 75,863 galaxies and represents both 
the sample with which we cross-match the CG sample and the basis 
of the pool from which we construct a control sample. The metal- 
licity sample is significantly smaller than our initial star forming 
sample because the metallicity calculation requires additional emis- 
sion line detections beyond those required to classify the galaxy 
as star forming. We take stellar mas^ values as calculated by the 
MPA/JHU group from fitting models to the 5-band SDSS DR7 pho- 
tometry 1^ These masses are, in general, in good agreement with 
spectroscopically derived values calculated by |Kauffinann et al.| 



^ Unless otherwise stated, all mass values in this work are stellar masses. 
^ A discussion of their method is available at: |http : //www . | 
[mpa-garching .mpg .de/SDSS/DRZ/mass.comp .htmll 



( |2003a| ). Star formation rates within the SDSS fibre (3") are taken 
from [Brinchmann et al.| ( |2b04| ), who use a set of 6 emission line 
fits to determine the SFR. Aperture corrected SFRs are also present 
in the catalogue. These SFRs are corrected from fibre to total val- 
ues by calculating the galaxy light not contained within the fibre; 
by fitting models to the photometry of the galaxy outside the fibre, 
the aperture correction can be eff'ectively made (See [Brinchmannj 
|et al.|2004 for a full discussion of their aperture correction method- 



ology).[Sa] 



ilimet al.| ( [2007l ) compare the results of this methodology 
with the SFRs obtained from UV flux, and find that for star forming 
galaxies the two methods agree very well, with no bias introduced 
by the aperture corrections. We use the aperture-corrected SFR val- 
ues for the rest of this work. 



2.4 Compact Group Sample 

Our parent CG sample consists of 3491 CG galaxies in 828 CGs 
(Mil). This sample is a refinement of the Catalogue A of M09. 
Typical intergalactic distances within the groups are a few tens of 
kpc, ranging up to 140 kpc. The M09 CG sample appUes a mod- 
ification of the original criteria laid out by |Hickson| ( |1982| ) to the 
SDSS DR6, with isolation, minimum number of galaxies, and den- 
sity requirements: 

(i) N{^m = 3) > 4 

(ii) > Wg 

(iii) iJLe < 26.0 mag arcsec"^ 

N(Am = 3) is the number of galaxies within 3 magnitudes of 
the brightest galaxy within the group, jUe is the surface brightness of 
the group, 6g is the angular size of the group, and 6^ is the size of 
the circle beyond which there are no galaxies within the 3 magni- 
tudes required for group membership. These criteria guarantee that 
there must be at least 4 galaxies in the group and the group itself 
must be separated by at least 3 times its own radius from any other 
equally bright galaxies. The surface brightness criterion guarantees 
the compactness of the group. Applying these criteria, M09 find 
2297 CGs, containing 9713 galaxies, down to a limiting magnitude 
of r = 18. However, there is no redshift criterion in the identifica- 
tion process for these galaxies; galaxies identified as a CG are done 
so based only on their photometry, regardless of whether or not they 
have concordant redshifts. Approximately 50% of the original 9713 
galaxies have reliable spectroscopic redshifts in the DR7 (Mil). 

In order to address the issue of false (due to projections) CGs 
in the M09 sample. Mil apply a statistical likelihood restriction 
to the master CG sample. Briefly, by looking at the probability 
density functions for the combinations of spectroscopic and pho- 
tometric redshifts available for CG galaxies found within a given 
group, a large fraction of interloping galaxies can be removed. If 
this process of interloper rejection results in the group failing the 
richness criterion used to define the sample, the entire group is re- 
jected. Mil estimate that there could be 20%-30% contamination 
remaining in the cleaned sample. Mil identified a distinction in the 
distribution of distances between the CG and the next nearest clus- 
ter or rich group. Approximately 50% of the cleaned CG sample 
can be classified as within < 1 Mpc of a rich group, using the 
|Tago et aT] ( |2010| ) catalogue compiled from the DR7 (see Figure 
5 in Mil). The other half of the sample is distributed at distances 
> 1 Mpc away from group structures. Placing a dividing line in 
the CGs at the minimum between the two distributions. Mil de- 
scribe two populations: 'embedded' and 'isolated' CGs. We adopt 
this terminology for the galaxies residing within those groups. Mil 
shows that the galaxies residing within the two density bins have 
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systematically different photometric properties (i.e., an increase in 
the fraction of blue galaxies within isolated CGs). 

We cross-match the galaxies belonging to either embedded or 
isolated CGs with our final metallicity and SFR sample to select 
only the galaxies residing in CGs that also have robust metallici- 
ties and SFRs. This cross-matching results in our final CG galaxy 
sample of 112 galaxies. Splitting the sample into embedded and 
isolated CGs results in 62 galaxies in isolated CGs, and 50 galax- 
ies in CGs embedded within a large scale structure. For a list of 
the embedded and isolated CG galaxies and their properties, see 
Tables [aT] & |A2| respectively. SDSS thumbnails of 4 random em- 
bedded CGs and 4 isolated CGs are presented in Figures [aT] & |A2| 
respectively. 



2.5 Matching to Controls 

In order to make an eff'ective comparison between our CG sample 
of galaxies and a control sample of galaxies, we must eliminate 
any major sources of bias. The three major sources of bias are stel- 
lar mass, redshift, and environment (Perez et al. 2009 ). As we are 
looking for an eff'ect that varies with environment, we cannot match 
in that property, but we wish to match a control sample to the CG 
sample in stellar mass and redshift, such that the distribution of 
stellar mass and redshift in the control match the distributions in 
the CG sample. Our control pool is defined as any galaxy in the 
final metallicity sample that is not flagged as belonging to a CG. 

Galaxies in our CG sample are matched to galaxies in the con- 
trol pool simultaneously in total stellar mass and redshift, in a sim- 
ilar way to |Ellison et al.| { [200 8b). Briefly, this algorithm finds the 
best simultaneous match in mass and redshift to each galaxy in the 
CG sample. Once every galaxy has been matched to a control, a 
Kolmogorov-Smirnov (KS) test is used on the total distributions of 
the CG sample and control. If the KS test finds that the distributions 
are consistent with being drawn from the same parent distribution 
at > 30%, then the matching continues without replacement until 
a maximum of 50 control galaxies have been matched to each CG 
galaxy in our sample, or the KS-test results in a probability < 30%. 
Since the pool of possible control galaxies is large compared to 
the CG galaxy sample, the matching procedure reaches the full ca- 
pacity of 50 matches per CG galaxy. Our final CG samples of 62 
isolated CGs and 50 embedded CGs thus have control samples of 
3100 and 2500 galaxies respectively. The masses and redshifts of 
the control galaxies are typically matched to within 0.05 dex and 
0.003 respectively of the CG galaxy value. The E(B-V) values of 
the galaxies are found to be consistent between all four samples 
when a KS-test is applied. 

The resulting normalised distributions of mass and redshift are 
shown in Figure [2] The control samples and the CG samples can 
be seen to match very well in both mass and redshift (the KS-test 
probabilities are > 99% between CG galaxies and their controls for 
both mass and redshift distributions). The embedded and isolated 
galaxies are also statistically similar to each other in both stellar 
mass and redshift, with KS-test values of 32.57% and 49.16% re- 
spectively. 



3 OFFSETS IN THE SFR AND METALLICITY OF 
COMPACT GROUP GALAXIES 

In this section, we present our methodology for finding and quan- 
tifying diff'erences between the SFRs and metallicities of compact 
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Figure 2. Normalized distributions of isolated (left panels) and embedded 
(right panels) CG galaxies and their controls for total stellar mass (upper 
panels) and redshift (lower panels). Galaxies are matched to controls in 
total stellar mass and redshift simultaneously until 50 control galaxies are 
matched to each group galaxy. KS-test probabilities that the samples are 
drawn from the same population are > 99% for the control and CG galaxy 
distributions in all panels. 



group galaxies and their controls. We also present a set of bootstrap 
simulations to quantify the significance of a given measured off'set. 



3.1 Offset methodology 

To quantify shifts away from the expected metallicities and SFRs 
at a given mass, we use a modification of the off'set method from 
[Patton et al.| ( |201l) . For each galaxy in the CG sample, we take 
the median (metallicity or SFR) value of all 50 control galaxies to 
which it was matched, and subtract that value from the observed 
value of the CG galaxy. An off'set is therefore defined as: 

Alog(0/H) = (log(0/H) + 12),,_,-^i/2 (log(0/H) + 12)_,,,,, , (4) 

Alog(SFR) = log(SFR),,,,,,,, - 1/2 log (SFR),,,,,,,,, (5) 

where yL/1/2 signifies the median. The calculated off'sets are pre- 
sented in Figures |3] and |4] The median metallicity ofi'set for the 
isolated CG galaxies (solid vertical line in Figure|3]) is -0.02+0.04 
dex (metal poor), whereas the embedded population has a median 
metallicity off'set of 0.00+0.02 dex (no off'set; dashed vertical line). 
The uncertainty on the median is calculated using a jackknife tech- 
nique. For each sample (embedded and isolated) the median off'set 
is re-calculated by systematically removing one galaxy. For a sam- 
ple of N galaxies, the jackknife error on the median is given by 



N 



(6) 



where /Ui/2(N) is the median value for the full sample of N galaxies, 
f^\/2(N - 1) is the median value fovN-l galaxies, and each of N 
galaxies is removed from the sample once. 

The SFR off'set distribution is shown in Figure |4j the embed- 
ded CG galaxies have SFRs that are lower than the control, with 
a median off'set of -0.03+0.05 dex. Isolated CG galaxies, by con- 
trast, have a median off'set of +0.07+0.03. Again, uncertainties are 
determined from the jackknife statistic. 
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Median^ -0.02± 0.04 

Median^ 0.0± 0.02 
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Metallicity offset (A log(0/H)) 

Figure 3. Offsets from the control MZR as defined by the median value per 
set of control galaxies for the embedded (red dashed line) and isolated (blue 
solid) compact group sample (50 and 62 galaxies, respectively). Black solid 
and dashed lines show the median offsets for the isolated and embedded 
compact group galaxies. Embedded galaxies show no median metallicity 
enrichment, contrasted with the median metallicity offset of isolated groups, 
which show a dilution of -0.02+0.04 dex compared to the overall relation. 



Median^ 0.07± 0.03 

Median^ -0.03± 0.05 

I I Isolated 

CZ2 Embedded 



SFR offset (A log(SFR)) 

Figure 4. Same as FigurejS] but now showing the offset of CG galaxies from 
the SFR - mass relation. Isolated CGs show a median SFR enhancement 
of +0.07±0.03 dex, while embedded CG galaxies have SFRs depressed by 
-0.03 ±0.05 dex below the control sample. 



In order to determine whether the measured median offsets 
are strongly influenced by the errors on the SFRs, we use a Monte- 
Carlo resampling simulation to test the sensitivity of our results. 
We define a distribution of possible SFR values as a gaussian with 
mean of the original SFR, and with width defined as the Icr er- 
rors on the SFRs, taken from B rinchmann et al.) ([2004 ). Median 
errors on any individual galaxy are ~ 0.1 dex. We draw new SFR 
values for both the CG and the control samples, and recalculate 
the offsets for the CG samples. We find that the median offsets 
are not significantly changed. After 10,000 iterations, the median 
A log(SFR) for isolated CG galaxies is 0.07 with a Icr scatter in 
the distribution of 0.05. Embedded CG galaxies show a median A 
log(SFR) of -0.06 + 0.06 (see Figure |A3]l. These simulated me- 
dian offsets and their uncertainties are consistent with the medians 
in Figure |4] and we conclude that the errors on the SFR measure- 
ments are not artificially enhancing the separation between the two 
samples. 

The uncertainties on the median offsets indicate that the metal- 
licities of the galaxies within both embedded and isolated CGs are 
consistent with a mass and redshift-matched control sample of non- 
CG galaxies. While the median SFR of the embedded CG galaxies 
is consistent with the control, the median SFR of the isolated CG 
galaxies is significantly offset from its control at ~ 2cr. Further, 
the isolated and embedded medians are significantly different from 
each other at the Icr level. 

3.2 Significance simulations for the SFRs 

The jackknife technique estimates the statistical error on our sam- 
ples of ~ 60 galaxies. However, we also want to quantify how of- 
ten a given median offset will occur by chance (due to the intrinsic 
scatter in the SFR-mass relation). A random sample of 60 galax- 
ies is therefore bootstrapped from the control sample 10,000 times. 
Recall that each control is one of 50 that has been matched to a 
given CG galaxy. When one of them is selected at random, the 49 



others in the set which are not selected as the 'sample' are left as 
the 'control'. Offsets are then calculated between the randomly se- 
lected galaxy and each of the remaining controls. The median offset 
for the sample of 60 galaxies is determined, along with the jack- 
knife error on that median. As the galaxies selected are part of the 
control pool, and no offset is artificially introduced, the majority of 
these offsets should be centred around zero. For each median, the 
significance cr is calculated by dividing the measured median off- 
set by the jackknife error on the median. The fraction of the total 
number of runs with significant offsets (defined as cr > 1, 2, 3, or 4) 
can then be calculated. 

We can now compute for any given offset the fraction of the 
total number of bootstraps which returned significant medians at or 
beyond the offset being investigated. The results are shown in the 
top panel of Figure [5] The different lines indicate different values 
of cr. To compare to our results, the SFR enhancement seen in iso- 
lated CGs would be along the 2cr line. This simulation indicates 
that there is a 0.19% chance of detecting a SFR enhancement of 
+0.07 +0.03 dex at random. We therefore conclude that the SFR 
enhancement we have detected, at its given observed significance, 
has a very low probability of occurring by chance due to scatter. 

Whilst the above procedure assesses the probability that a 
given offset is found between the CG and control samples, we can 
also quantify the significance of the difference between the embed- 
ded and isolated samples. This is achieved by repeating the boot- 
strapping procedure described above, but this time selecting 2 sep- 
arate samples of 60, and the absolute value of the difference in their 
median offsets is found. We then find the percentage of iterations 
where the simulated differences in median is greater than a given 
offset. The results of the two sample bootstraps are shown in the 
lower panel of Figure [5] The probability of a range of 0.10 dex be- 
tween two samples drawn at random from the control SFR-mass 
relation is 0.67%, which corresponds to > 2cr significance. 
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Figure 6. Results of the sensitivity simulation for the MZR offsets. Each 
point is the result of an independent simulation, bootstrapped 5,000 times. 
Each simulation selects a random sample of 60 galaxies from the control 
sample, and inserts an artificial off'set i into the metallicities. The off'sets are 
then calculated as normal, resulting in a measured off'set m. The fraction of 
those runs which result in significant detections at the 1,2, and 3 cr level are 
plotted in black, red, and blue, respectively. The median diff'erence between 
the measured off^set m and the input off'set / for the fraction of significant 
runs (at 3cr) is displayed as the blue squares in the lower panel as a function 
of i. Jackknife errors on the medians are smaller than the points. 



Figure 5. Results of the bootstrapping simulation for SFRs. A randomly 
selected sample of 60 galaxies is pulled from the control pool, and their 
off'sets from their matched controls recorded, repeated 10,000 times. The 
top panel is the fraction of the total number of runs where a significant off'set 
(i.e., cr > 1, 2, 3, or 4) greater than a given off'set i is found. Off'sets of -0.04 
and +0.08 dex at 2 and 4cr significance have 2.14% and 0.02 % likelihood 
of arising by chance, respectively. The bottom panel tests the percentage 
of times that a range in medians between samplings occurs; a diff'erence in 
median of 0.10 dex could occur randomly ~ 0.67% of the time (> 2cr). 

3.3 Upper limit simulations for metallicities 

Although we find a statistically significant diflference between the 
SFRs of embedded and isolated CGs, the metallicity oflfsets are in- 
conclusive (-0.02 + 0.04 and 0.00 + 0.02 for isolated and embedded 
CGs respectively). Metallicity oflfsets in both embedded and iso- 
lated structures are statistically consistent with their control sam- 
ples. To test the sensitivity of our null result given the sample size, 
we run a set of simulations similar to the bootstrap used in ^3.2| but 
now introducing an artificial oflfset. In this simulation, 60 galaxies 
are randomly selected from the control sample and an artificial, 
single value oflfset (/) is added to all the metallicities of the selected 
galaxies. The resultant metallicities are then run through the oflfset 
calculation algorithm. Oflfsets are calculated between each of the 
60 randomly selected galaxies and the remaining 49 in its matched 
set. The median measured oflfset (m) of the 60 galaxies is recorded. 
This step allows us to quantify the diflference between the input oflf- 
set / and the recovered oflfset m through our algorithm. We expect 
/ - m = if the sample can accurately recover the input metallicity 
oflfset. This bootstrap process is repeated 5,000 times. We calcu- 



late the jackknife error on each of the median values returned by 
the bootstrap. For each realisation, we calculate a significance cr 
by dividing the measured oflfset m by the jackknife error. We can 
then find the fraction of the total number of bootstrapped samples 
which return significant oflfsets for a given /. This fraction is then 
recorded, and plotted in the top panel of Figure [6] Red, green, and 
blue points reflect the fractions for cr > 1,2, and 3, respectively. 
Each inserted oflfset requires a separate, independent run of 5,000 
iterations of the bootstrap, and is plotted as a separate point in the 
figure. 

In the bottom panel of Figure[6] the median diflference between 
the input and measured oflfsets (/-m) is shown for those runs which 
returned significant medians at the 3cr level, as a function of the in- 
serted oflfset /. This panel indicates that we are able to accurately 
recover an oflFset, if it is systematically diflferent from the control. 
For values of / smaller than +0.05, (/-m) shows an asymmetric shift 
away from zero. At very small /, it becomes increasingly unlikely 
for an oflfset m to be measured as significant. Therefore, those sim- 
ulations which are measured as significant at small / are likely to 
be those which have scattered to median values further from zero. 
Negative values of / therefore preferentially return more negative 
values of m, and positive values of / return more positive values of 
m. 

With these simulations as a base, we may compare the frac- 
tions of significant measured oflfsets returned by any given iteration 
of the bootstrap with the magnitude of oflfset we might expect to 
see. We use what has been found in previous studies of metallicity 
changes in high density environments as a guideline for reasonable 
magnitudes of oflfsets. The metallicity suppressions found in the 
pairs sample of [Ellison et al.| p008b| ) and of |Michel-Dansac et al.| 



© 0000 RAS, MNRAS 000, 000-000 



8 Scudder et al. 



( |2008| ) are of order 0.05 dex, and the enhancements in overdensi- 
ties in [Cooper et al.|([200 8 ) are ~0.03 dex. Similar magnitudes are 
found by lElHson et al.| ( |2009j ), who find off'sets of 0.04 dex in galax- 
ies in high density regions, relative to their control. [Mouhcine et al.| 
P007^ find slightly higher diff'erences in metallicity of 0.06-0.08 
dex between high and low density environments. For off'sets in the 
range 0.03 - 0.05 dex, we find we only have a 34.8-51.1% chance 
of detecting a 3cr result. We do not reach 99% confidence levels of 
detecting a Icr off'set with our sample size until we reach off'sets of 
0.13 dex, which is much larger than anything previously discussed 
in the literature. The lack of significant metallicity off'set in our 
CG sample may therefore be a limitation of our sample size. The 
fact that our metallicity results are both consistent with the control, 
combined with the results of this simulation, indicate our sample is 
unlikely to result in a significant off'set of the magnitude suggested 
by the literature. 

In summary, the off'sets in SFR in our CG samples appear to be 
statistically distinct from each other, and galaxies in isolated CGs 
have significantly enhanced SFRs, whereas embedded CGs show 
no enhancement. However, the measured metallicity off'sets are not 
sufficiently sensitive to probe metallicity differences between em- 
bedded and isolated compact group galaxies and their controls. 



4 DISCUSSION 

We have used the varying large scale environments of a consistently 
selected sample of compact, high density groups to probe the influ- 
ence of the presence of a rich group or cluster. The compact group 
selection criteria allow us to select a consistently dense local en- 
vironment on scales of ~ 100 kpc (see Figure 6 of Mil) which, 
ignoring all external influences, ought to be subject to the same 
internal evolutionary processes. Galaxies in CG systems typically 
have small projected separations from each other; for our sample 
the median projected separation to the nearest neighbour is 46 kpc 
for both isolated and embedded samples. In galaxy pairs separa- 
tions < 30 kpc result in significantly enhanced SFRs (e.g.. Barton^ 
etal . 2000, Lambas et al. 2003, Alon so et al.|2004||2506l|Emson 
et al.|2008bn2010t|Scudder et al. 2012^. Our sample has a signif- 
icant population of galaxies with separations of this order, making 
galaxy-galaxy interactions an appealing explanation for the SFR 
enhancement seen in the isolated CG sample. The level of offset 
seen in our sample of isolated CG galaxies (+0.07 dex) may appear 
modest, but is comparable to studies of merger-induced star for- 
mation, which generally find median SFR enhancements of around 
+0.1 dex (e.g., |Emson et al.|2008b| [Scudder et ari20r2 ). If the iso- 
lation criterion of the CG selection is eff'ective at limiting outside 
influences, and interactions are the only dominant mechanism for 
transformations, we might expect both samples of CG galaxies to 
show consistent SFR enhancements over the control. 

We find that galaxies in CGs in close proximity to a cluster 
or rich group (<1 Mpc) do not show the same preferential en- 
hancement as their counterparts in isolated CGs. The embedded CG 
galaxies instead show consistency with the average galaxy in the 
control. Galaxies in isolated CGs (>1 Mpc from large scale struc- 
tures), by contrast, show SFR enhancements over the control, con- 
sistent with the general trend presented by the SFR-density relation. 
This SFR disparity indicates that the local overdensities identified 
by the CG selection criteria cannot be considered in isolation. The 
presence of a nearby rich group exerts some influence on the star 
forming galaxies found within CGs, either by introducing an ad- 
ditional quenching mechanism, or because the galaxies within rich 



groups are physically distinct from their isolated counterparts. The 
SFRs and metallicities of the galaxies within those groups therefore 
serve as diagnostics of the impact of the larger scale environment. 
However, given the small magnitude of the difference between the 
two CG samples, the large scale environment does not appear to 
drive major changes in a galaxy's evolution at fixed local density. 

4.1 Density scale dependences 

Having determined the sensitivity and significance of the metallic- 
ity and SFR off'sets, we wish to compare these results to previous 
work. However, previous metrics of environment usually fall into 
one of two categories; either an n^^ nearest neighbour approach, 
or calculating galaxy density within annuli of increasing radius 
around a given point. In this work, we have used a new method 
of determining the large scale environment of a galaxy. Our divi- 
sion of galaxies into embedded and isolated systems is a binary as- 
signment, based upon the presence (or lack) of rich group structure 
within a cylinder in redshift and projected separation space, centred 
upon the CG. This is an entirely diff'erent method of characterising 
the nature of a galaxy's environment than has been typically used 
in the past. 

Comparisons to the results found using existing measures of 
environment will require attention to the methodology of those 
metrics to be sure that they can probe the same kinds of overden- 
sities. Since the n^^ neighbour and annulus methods are calculated 
using entirely different algorithms, the meaning of the densities re- 
turned (as it applies to our CG sample) will vary, and comparisons 
to the results presented here must bear this in mind. In addition, 
most previous work has focused upon the relative importance of 
local and large scale dependences. This relative comparison is one 
that we cannot address in this work, as our classifications of em- 
bedded and isolated are not relative selections, and the local den- 
sity - that of the CG itself - is relatively constant across our sample. 
Therefore, in comparing to previous work, we can only determine 
whether or not the previously obtained results are consistent with 
some level of SFR dependence on large scale environment as mea- 
sured by our own method, to the extent that the density metrics may 
be compared. 

The n^^^ nearest neighbour approach identifies high density re- 
gions of space by returning a shorter distance to the n^^ neigh- 
bour in dense regions, and a longer distance in lower density re- 
gions (e.g., [Hashimoto et al. |1 998 1 [Carter et al.|2001| ). Varying the 
value of n eff'ectively changes the distance scales over which one 
wishes to probe the densities. With small values of n, the densi- 
ties returned will be more reflective of local scale densities. Larger 
values of n will necessarily wander further afield to reach more 
distant neighbours. These distant neighbours will return a much 
larger volume of space over which to average the galaxy density, 
relative to their locally dense counterparts. The larger distances for 
isolated systems means that the densities derived from the iden- 
tified n^^ neighbour will be more reflective of a large scale den- 
sity than those in locally dense environments (which return smaller 
distances), even if those local overdensities are relatively isolated 
from other structures ( [Weinmann et al.[|20 06). Since the distance 
to n^^ nearest neighbour can vary smoothly from system to system, 
the distance scales being probed will also vary strongly as a func- 
tion of local galaxy density. Furthermore, since the distance to n^^ 
nearest neighbour is converted into a density, information about the 
distribution of those galaxies on scales less than that distance is 
lost. For the purposes of comparison to our metric, the n^^ nearest 
neighbour algorithm simply needs to be able to distinguish between 
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isolated and embedded CG systems. Both isolated and embedded 
CGs have similar densities in the cores, but exhibit different galaxy 
surface densities on scales of 0.1 - 1 Mpc (Mil). As long as the 
distances returned by the n^'^ nearest neighbour probe these scales, 
they should be able to distinguish between the isolated and embed- 
ded systems. 

This relatively straightforward criterion for comparison with 
the n^'^ nearest neighbour provides a sharp contrast with the other 
main metric of galaxy density. The annulus approach is more diffi- 
cult to compare to our method. In this approach, the galaxy density 
within annuli of increasing radius is measured (e.g., [Lewis et al.| 
[20021 [Kauffmann eraLl[2004l [Blanton et aLl[2056t . Many works 
place their dividing line between 'local' and 'large' scale at 1 Mpc. 
The densities within the annulus of radius 1-3 or 1-6 Mpc is there- 
fore used as their measure of large scale environment. 1 Mpc is 
also a natural dividing line in our sample, but only in the sense 
that this is the distance to the nearest group that distinguishes iso- 
lated and embedded CGs. The annulus approach is sensitive to the 
presence of broad overdensities beyond the distance of the defined 
'large scale', generally 1-6 Mpc away from the selected galaxy. 
Our method provides the ability to define a system as embedded 
even when the overdensity is entirely on scales < 1 Mpc. In this 
case, the CG would still distinctly be embedded in a more massive 
halo than its own, but would not be distinguished as 'large scale' 
structure by an annulus of radius >1 Mpc. If some fraction of the 
rich group extends beyond 1 Mpc distant from the CG galaxy, then 
it would be detected by the >1 Mpc annulus. However, Mil shows 
that rich groups and CGs in any environment converge to similar 
galaxy densities at separations ~ 1 Mpc. With the median distance 
from CG centre to a rich group for embedded groups being of order 
~ 0.3 Mpc, and the majority of rich groups 's strong galaxy surface 
density visible on scales < 1 Mpc, we might expect that only a small 
fraction of the galaxies within a rich group to extend beyond the 1 
Mpc threshold. Furthermore, the annulus method loses all infor- 
mation on scales smaller than the area of the annulus. If an annulus 
contains an overdensity, but also contains a very low density region, 
these two extremes will average to a value that is not representative 
of either region. 

With these differences in mind, we can now attempt to com- 
pare our results of a SFR dependence on large scale environment 
with previous studies. We emphasize that we are investigating only 
the dependence on those galaxies which are star forming relative to 
a star forming control, not the overall trend of star formation rate 
with density. Those works which conclude that the large scale en- 
vironment affects the evolution of star forming galaxies by way of 
the cluster-centric radius ( Park & Hwang 2009 ) or dark matter halo 
mass ( [Goto et al.[2003[[Weinmann et al. 2006j are compatible with 
our results. However, a large body of work in the literature (using 
a mixture of methodologies) has concluded that the local scale ef- 
fects dominate galaxy evolution (e.g., Hashimoto et al.'1998 , Carter 
et al. 2001 , Lewis et al. 2002 , Kauffmann et al. 2004 , Blant on et aTj 
|2006^ Blanton & Berlind^2007j . These analyses consistently find 
a strong SFR dependence upon local galaxy density, with either 
very weak or absent dependence upon galaxy densities beyond the 
1 Mpc scale. For example, Blanton et al. ( 2006) constrains the ef- 
fects of large scale (1-6 Mpc) effects to be small, but does not ex- 
clude them as a minor contributor. Our results are not in disagree- 
ment with those studies which found that the local scales dominate 
galaxy evolution, consigning the large scale effects to a minor role. 
The small difference in A log(SFR) between star forming galaxies 
in isolated and embedded CGs surely confirms the minor nature of 
the large scale environment. Our work indicates simply that there is 



a non-zero effect upon the galaxies within a larger structure beyond 
local density effects; it does not require that the effect must be the 
dominant one. Other works in the literature are less conservative, 
and interpret the lack of a strong dependence upon the galaxy den- 
sities at > 1 Mpc radii as evidence that the large scale environment 
does not have an effect upon a galaxy. Blanto n^ Berlind| ( [2007| ) 
find that for the star forming fraction, large scale structure (defined 
as structure in an annulus > 1 Mpc) has no discernable effect upon 
the galaxy, with a bias in their measurement of no more than 5%, 
and density error of order 15%. Kauffmann et al. ( 2004) found that 
local densities dominate the SFR evolution; at low masses, the SFR 
of a galaxy was up to 10 times lower in high local densities than it 
was in low local densities. When testing the spatial scales between 
1-3 Mpc, they found a very noisy relation in the 4000 A break 
strength, their metric for SFR, with density: consistent with no de- 
pendence. In contrast, we have found that SFRs do depend, albeit 
at a low level, on their association with a large scale structure. 

Reconciling our results with those in the literature likely 
comes down to the nature of the 'environment' identified by the 
different density metrics that have been previously used. Our divi- 
sion in environment identifies galaxies which are within 1 Mpc of 
a larger group than the CG itself. By the time scales of > 1 Mpc 
from the centre of a group are reached, the galaxy densities for iso- 
lated, embedded, and rich groups have all converged (Mil). This 
indicates that most of the galaxy overdensity due to the rich group 
has also dwindled by 1 Mpc from the centre. If the groups are then 
sufficiently far away from the CG, they might cross over into the > 
1 Mpc annulus, and be detected as an overdensity by that metric. 
However, the peak of the distribution of the distance to the near- 
est rich group is ~0.3 Mpc for embedded structures (Mil). At 0.3 
Mpc, the majority of the rich group will be found within 1 Mpc of 
the CG, and would not significantly increase the measured density 
for annuli >1 Mpc. Therefore, even though the embedded CGs are 
situated in a significantly different environment than the isolated 
CGs, the overall galaxy density no longer reflects this difference at 
scales > 1 Mpc. Dividing environmental factors into large and local 
scale effects simply by using a cut in radius may blend extremely 
local drivers of evolution and some larger scale environmental fac- 
tors. Our metric therefore cleanly separates those groups embedded 
in a more massive halo from those which are independent halos, no 
matter what the projected scale of the more massive group may be. 

4.2 Physical Drivers 

The interpretation of the difference between embedded and isolated 
galaxy SFRs depends on whether the galaxies are intrinsically and 
physically different when they are found within an overdensity, rel- 
ative to in isolation, or whether the overdensity is imposing an addi- 
tional processing mechanism on an initially similar set of galaxies. 
If the former, then galaxies found within large scale structure are 
reacting to the same set of evolutionary processes as galaxies in 
isolation, but their intrinsic properties make SFR enhancement (for 
example) more difficult to achieve. 

One mechanism by which a physical difference in the galaxy 
populations between embedded and isolated CGs might arise is as- 
sembly bias (e.g., [Croton et al.|2007[ [Cooper et al.|2 010). Assem- 
bly bias states that galaxies in high mass potentials tend to be older, 
as they collapsed and began to cluster at earlier times. In this sce- 
nario, one might expect clusters to contain a higher fraction of old, 
so-called 'red and dead' galaxies, galaxies which have undergone 
many mergers, or galaxies which have been affected by the clus- 
ter potential through quenching mechanisms for a longer period 
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of time. Galaxies in rich groups tend to be preferentially gas poor 
(e.g., [Chung et al.|2009| as a result. This framework offers one ex- 
planation for the SFR suppression in overdense regions. If galaxies 
within rich groups are simply older than their counterparts in the 
isolated groups, then they would be further along in their evolu- 
tion, and may be further along the process of gas exhaustion and 
consumption. If the CG galaxies, due to their age, in addition to 
any environmental effects, are preferentially gas poor, then it will 
be significantly more difficult to trigger a SFR enhancement af- 
ter an interaction (e.g., Alonso et al.|2006| . Their isolated counter- 
parts, having collapsed at later times, are younger, more gas-rich, 
and at a much earlier stage in their evolution; these galaxies will 
show SFR enhancements much more easily, as their gas reservoir 
is less depleted. This explanation does not require a difference in 
the mechanisms that modulate star formation rates in the two en- 
vironments; the difference in intrinsic galaxy properties is enough 
to alter the galaxy's response to the same large scale environmental 
mechanisms. However, all the galaxies in our sample are emission 
line galaxies. While their star formation rates may be lower than 
the average galaxy, they still have significant amounts of gas. They 
cannot have been so affected by quenching and stripping mecha- 
nisms as to remove all their gas, if they can maintain strong emis- 
sion lines. A large portion of the galaxies in our sample, regard- 
less of their classification as embedded or isolated, retain relatively 
late-type morphologies (see Figures [aT] & |A2| ). The maintenance 
of the emission lines and their late-type morphologies limits the 
magnitude of the environmental and age-related evolution that the 
galaxies could have undergone. 

Turning to consider processes that could alter the SFRs of a 
given galaxy in a different environment, we can eliminate quench- 
ing mechanisms that function on very small scales such as those 
associated with interactions; these should be consistently affect- 
ing both isolated and embedded systems. Whatever mechanism is 
preventing SFR enhancement in embedded CG galaxies must be 
associated with the larger scale cluster potential, and not simply 
with the high densities of the CG system. One explanation could 
be the difference in evolutionary processes in central and satel- 
lite galaxies. Galaxy populations divided in this way show that 
galaxies which are not the dominant galaxy in the dark matter halo 
undergo stronger environmental processes than the central galaxy 
does ([Haines et al.|2006|[Kimm et al.|2009||Peng et al.|2011[rnn^ 
[ker et al.|2011| ). Embedded CGs would, by definition, be a satellite 
halo within the larger cluster halo. The enhanced processing seen 
in most satellites would therefore be systematically applied to all 
galaxies in the embedded systems. [Blanton & Berlind[ ( [2007| ) find 
that high density substructures (on scales of < 300 kpc) within a 
cluster are more correlated than expected with hosting red galax- 
ies, i.e., these galaxies are further along in the quenching process. 
By contrast, the probability of the dominant galaxy in the poten- 
tial being in our sample is significantly higher for isolated systems, 
where the CG environment is more likely to be an isolated halo. 

If the difference between galaxies in isolated and embedded 
CG systems is not due to a systematic physical difference in galaxy 
population, then the galaxies in both systems ought to have begun 
their evolution in the CG as a roughly uniform population of galax- 
ies. In order to explain the absence of the SFR enhancement seen 
in the isolated CG galaxies, an additional evolutionary force must 
be at work in embedded CG galaxies due to the cluster potential 
that is not in effect in isolated halos. This brings into suspicion 
the main drivers of galaxy evolution in clusters: ram-pressure strip- 
ping, galaxy-galaxy harassment, and strangulation. The isolation 
criterion for CG selection means that external galaxy harassment 



is unlikely to be a significant influence. Galaxy harassment within 
the CG would contribute similar levels of evolutionary processing 
to both galaxy populations of CGs, and is therefore not a useful 
explanation for the difference in SFR response. Further, CGs have 
relatively low velocity differences, so harassment within the group 
itself is unlikely to dominate over lower velocity tidal interactions. 
This leaves ram-pressure stripping (removal of cold gas from the 
disk) and strangulation (removal of the external gas supply). [Wein-[ 
[mann et aL] ( [2006| ) argue that ram-pressure stripping should af- 
fect low mass galaxies more quickly, as they have a smaller po- 
tential and are more easily stripped of their gas, within a halo of 
given mass. This ought to result in a higher fraction of low mass 
early type galaxies relative to the fraction of high mass early types. 
Their sample, taken from the SDSS DR2, contradicts this predic- 
tion, finding no evidence for preferential changes to the lower mass 
galaxies. Wetzel et al.[ ( [201 1) find that galaxies introduced into a 
larger halo can retain their SFR at pre-accretion levels for 1-2 Gyr 
before rapid quenching takes place; if ram-pressure stripping were 
the dominant mechanism, the galaxies would have their cold gas 
removed on much quicker timescales. Ram-pressure stripping is 
also excluded by a number of other authors (e.g., jTanaka et al.[ 
[2004} [Blanton & BerHnd[[2007] ) due to the large radial extent of 
SFR suppressions, and SFR suppressions in low mass halos. Nei- 
ther of these environments are expected to have the dense intra- 
cluster medium needed to support effective ram-pressure stripping 
(e.g., Rasmussen et al.|2008j ). This leaves strangulation as the prime 
suspect for SFR suppression in cluster environments, as it could 
still remove the hot halo of gas around a galaxy in lower density 
or lower mass environments (Balogh et al. 2000 ). Furthermore, it 
does not affect the morphologies or cold gas reservoirs of the galax- 
ies involved ( [Weinmann et al.|2006| ). Mil find that the morpholo- 
gies within embedded and isolated CGs are very similar, which also 
supports the idea of an environmental process being at work, rather 
than different evolutions due to age or gas exhaustion. Strangula- 
tion therefore appears to be an appealing mechanism for long-term 
suppression of SFR in emission line galaxies. 



5 CONCLUSIONS 

We present a sample of robustly calculated gas-phase metallicities 
and star formation rates for a sample of 75,863 star-forming galax- 
ies in the SDSS DR7. These metallicities use emission line fluxes 
which passed rigorous quality control, and are cleaned of any AGN 
contribution. From this sample, we define a compact group sample 
of 1 12 galaxies that is uniformly selected on criteria for local rich- 
ness and cleaned of likely interlopers. The compact group sample is 
further split by large scale environment into isolated and embedded 
subsamples, with 62 and 50 galaxies in each sample respectively. 
These two samples of star forming CG galaxies therefore have the 
same (very high) local density, but different large scale environ- 
ments. The samples of compact group galaxies are simultaneously 
matched in mass and redshift with 50 non-compact group galax- 
ies each in order to form robust control samples. By comparing the 
control samples with the compact group galaxy properties, we can 
measure the difference in the SFR and metallicity at fixed mass. 
The main conclusions of our work are as follows: 

(i) Isolated and embedded compact group galaxies are offset 
from each other, and isolated CG galaxies are significantly en- 
hanced relative to the field SFR-mass relation. Embedded galaxies 
do not show a similar enhancement. The SFRs of the galaxies in 
isolated CG systems show a median enhancement of +0.07+0.03 
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dex, whereas the galaxies in embedded systems show a median off- 
set of -0.03+0.05 dex. A series of bootstrap and Monte Carlo sim- 
ulations indicate that this offset is robust. Our results indicate that, 
even though previous works have found a primary dependence on 
local density, the SFRs of star forming galaxies in locally overdense 
regions are mildly sensitive to large scale environment at fixed local 
density. 

(ii) We find no evidence for a significant metallicity offset in 
either the isolated or embedded samples relative to their control 
samples. However, this non-detection is likely to be a result of small 
number statistics, as our sample is not likely to detect offsets of 
order 0.03-0.05 dex, which is the typical metallicity offset expected 
due to interactions or density driven effects. Our sample is therefore 
inconclusive with regards to the metallicity offsets. The simulations 
indicate that at 99% confidence we could have detected an offset of 
0.13 dex. 

Criteria classically used to select CGs define a consistent set of 
high density groups of galaxies. However, when these systems are 
divided according to the presence of a nearby cluster or rich group, 
the SFRs of the star forming galaxies within these systems become 
divided. The environmental split in median SFR demonstrates that 
large scale structure can exert an effect upon substructure found 
within it, even if at a low level, and that substructure is relatively 
isolated from other galaxies. The two possible explanations are: 
galaxies found within overdense regions are intrinsically different 
from galaxies in low density regions, which alters their response 
to the same local overdensities, or that the rich group structure im- 
poses an additional evolutionary mechanism onto the galaxies in 
embedded CGs. This may give further support to the differences 
in satellite and central galaxy evolution, with strangulation, found 
only in cluster environments, as an effective process in slowly shut- 
ting down star formation. 
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APPENDIX A: FIGURES & TABLES 

In the Appendix, we include additional figures and tables for illus- 
trative purposes. Figures [aT] and [A2| show Sloan images of sam- 
ple CGs from the embedded and isolated samples respectively. Ta- 
bles |aT] and |A2] contain information on the embedded and isolated 
CG galaxies respectively. Figur e] A3 [ shows the results of the Monte 
Carlo simulation conducted in |3| 
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Figure Al. 4 embedded CG galaxies, selected at random. The images are centred on the galaxy in the spectroscopic sample (identified by SDSS objID in the 
lower left corner), and are approximately 200 arc seconds to a side. 
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Figure A2. 4 isolated CG galaxies, selected at random. The images are centred on the galaxy contained in the spectroscopic sample (identified by SDSS objID 
in the lower left corner), and are approximately 200 arc seconds to a side. 
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SDSSObjID CGObjID Redshift Stellar Mass Metallicity SFR 



z log(M,) (Mo) 12+log(0/H) log(SFR)(Moyr-i) 
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0.072 


10.09 


9.09 
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9.64 


9.11 


-0.22 


587735348025819263 


SDSSCGAOl 108.4 


0.054 


9.81 


9.15 


0.06 



Table Al. Embedded Compact Group galaxy properties. SDSS Objid is the unique identifier for that galaxy within Sloan, CG Objid is the Galaxy ID from 
the M09 catalogue where numbers before the decimal indicate the group ID, and the number after the decimal varies between galaxies within one group. 
Metallicities calculations are described in %23\ SFR values come from the calculations of |Brinchmann et al.] ( |2004) . 
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SDSSObjID CGObjID Redshift Stellar Mass Metallicity SFR 

z log(M,) (Mo) 12+log(0/H) log(SFR)(Moyr-i) 



587739844854939883 


SDSSCGAOl 177.2 


0.114 


10.61 


9.10 


0.82 


588023720707424428 


SDSSCGA01088.4 


0.061 


10.11 


8.98 


-0.51 


587736975813836873 


SDSSCGA00878.2 


0.030 


9.97 


9.22 


0.72 


587739405707116711 


SDSSCGAOl 128.1 


0.054 


10.48 


9.04 


0.41 


587732471457186108 


SDSSCGA00146.1 


0.041 


10.24 


9.21 


0.13 


587725552264479032 


SDSSCGA01533.2 


0.066 


10.16 


9.18 


0.35 


587726101490827509 


SDSSCGA00370.4 


0.085 


10.25 


9.23 


0.57 


587725980687925438 


SDSSCGA02231.4 


0.097 


10.27 


9.18 


0.50 


588017724946972875 


SDSSCGA00942.4 


0.083 


9.79 


8.87 


-0.07 


587735344792928602 


SDSSCGA00422.5 


0.046 


9.23 


8.92 


-0.48 


587732702864080970 


SDSSCGA02097.1 


0.052 


9.67 


8.87 


0.39 


587736543623446530 


SDSSCGA00970.4 


0.087 


9.89 


8.97 


0.57 


587725041698799708 


SDSSCGA01076.4 


0.063 


9.61 


8.87 


0.12 


587725041698799706 


SDSSCGA01076.3 


0.064 


9.62 


8.87 


0.18 


587729386610557106 


SDSSCGAOl 836.3 


0.072 


10.08 


9.09 


-0.41 


587739377773576383 


SDSSCGA00513.3 


0.079 


10.13 


9.19 


0.49 


587728930803744828 


SDSSCGA00362.1 


0.112 


10.05 


9.06 


0.74 


587745244704473323 


SDSSCGA01087.3 


0.045 


9.84 


9.03 


-0.19 


587742062678900857 


SDSSCGA01309.3 


0.068 


9.69 


8.80 


0.07 


587745245240295524 


SDSSCGA00601.3 


0.113 


10.36 


9.20 


0.48 


587739827674022085 


SDSSCGA00041.2 


0.072 


10.39 


9.23 


0.82 


587729386068705391 


SDSSCGA00078.3 


0.056 


9.37 


8.84 


0.01 


587739844310007987 


SDSSCGA01088.3 


0.061 


9.93 


9.02 


0.06 


587731886268088671 


SDSSCGA01783.2 


0.078 


10.75 


9.25 


1.16 


587745539432776118 


SDSSCGA02092.1 


0.052 


9.48 


8.80 


-0.11 


587742774022045799 


SDSSCGA01953.4 


0.104 


9.88 


8.99 


0.34 


587742550686368002 


SDSSCGA00687.2 


0.087 


10.18 


9.19 


0.34 


588017627221721269 


SDSSCGA01937.2 


0.074 


10.25 


9.07 


0.64 


588848899373072598 


SDSSCGA00713.3 


0.118 


10.50 


9.01 


0.50 


588017990162710837 


SDSSCGA00718.2 


0.139 


10.84 


9.32 


1.11 


587732483824746731 


SDSSCGA00305.2 


0.084 


9.76 


8.88 


0.25 


587742013820305436 


SDSSCGA01570.2 


0.106 


9.95 


8.90 


0.84 


587742577521000590 


SDSSCGAOl 809.3 


0.093 


10.33 


9.19 


0.48 


588017627221721251 


SDSSCGA01937.3 


0.074 


10.09 


9.11 


0.53 


588007006048485395 


SDSSCGA00936.1 


0.074 


10.19 


9.07 


0.80 


587736975813836874 


SDSSCGA00878.3 


0.031 


9.71 


9.16 


-0.50 


587738947198910752 


SDSSCGAOl 170.2 


0.043 


9.71 


8.94 


-0.19 


587734893827129437 


SDSSCGA00945.1 


0.060 


10.56 


9.16 


0.79 


587741720678826052 


SDSSCGA00562.2 


0.067 


9.50 


8.90 


0.39 


587726101490827456 


SDSSCGA00370.3 


0.086 


10.20 


9.00 


0.56 


587725469057351815 


SDSSCGA01097.1 


0.132 


11.09 


8.89 


0.87 


587736782540964307 


SDSSCGA01542.3 


0.124 


10.69 


9.27 


0.76 


587739115237015957 


SDSSCGA01970.3 


0.099 


9.74 


8.88 


0.21 


587739610230554751 


SDSSCGA00820.2 


0.073 


10.20 


9.09 


0.37 


588017627221721162 


SDSSCGA01937.4 


0.074 


10.40 


9.03 


1.03 


588848900472111426 


SDSSCGA01473.2 


0.076 


10.11 


9.18 


0.48 


587739405707116699 


SDSSCGAOl 128.2 


0.054 


10.46 


9.08 


0.67 


588848899392340185 


SDSSCGA00476.3 


0.093 


10.44 


9.21 


0.62 


587729407546228763 


SDSSCGA02056.3 


0.029 


9.82 


9.02 


-0.46 


587739811028992231 


SDSSCGA01347.2 


0.068 


9.57 


8.93 


0.23 


587722984441119019 


SDSSCGA01703.2 


0.120 


10.36 


9.13 


0.74 


587722984441119010 


SDSSCGA01703.3 


0.121 


10.68 


9.27 


0.72 


587736546852864070 


SDSSCGA00337.4 


0.069 


9.79 


9.01 


0.49 


587737810644697402 


SDSSCGA01583.3 


0.103 


9.81 


8.95 


0.46 


587742578054201471 


SDSSCGA01218.3 


0.099 


10.53 


9.20 


0.52 


588016878825701672 


SDSSCGA00227.3 


0.082 


10.12 


9.05 


0.71 


587739610230554752 


SDSSCGA00820.3 


0.072 


9.91 


8.92 


0.35 


588009368008261825 


SDSSCGA02225.5 


0.092 


10.13 


9.04 


0.44 


587738411408294053 


SDSSCGA00614.3 


0.079 


10.00 


8.93 


0.33 


587742013297131696 


SDSSCGA01348.3 


0.102 


9.96 


8.96 


0.54 


587741603095707807 


SDSSCGA00610.3 


0.081 


10.08 


8.99 


0.19 


587742774555508821 


SDSSCGAOl 132.3 


0.123 


10.22 


9.14 


0.59 



Table A2. Same as Table |Al| but for the isolated CG galaxies 
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